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1 Introduction 

Nonlinear evolution equations in mathematical physics are largely stud- 
ied because of their important roles in most of the branches of science. The 
Klein-Gordon (KG) equation is a relativistic version of the Schrodinger's 
equation which describes field equation for scalar particles (spin-0). The 
KG equation has been the most frequently studied equation for describing 
the particle dynamics in quantum field theory. This equation with various 
types of potentials has appeared in field theories which is called nonlin- 
ear Klein-Gordon (NKG) model. Pion form factor [Ij, Hadronic atoms |2], 
Josephson junction [3j, Superfluid current disruption [4j, Shallow water [5] 
and Hawking radiation from a black hole |6J are some examples of the NKG 
applications. 

In recent years several solitonic solutions for the NKG equations have 
been proposed using different methods. These solutions have appeared in 
homogeneous and well behaved medium while in real world the medium 
of propagation contains disorders and impurities which add local space- 
dependent potentials to the problem. The behaviour of the NKG solitons 
during the interaction with these local potentials and their stability after 
the interaction are important questions. Therefore, investigating these sit- 
uations is very interesting. It is a very important subject because of its 
applications and also due to mathematical point of view. Motivated by 
this situation, we have studied the interaction of the NKG solitons with 
defects using different methods and the results are presented in this paper. 
Therefore a brief description of the NKG model and its solitons has been 
introduced in section 2. The motivation of the present work and some ap- 
plications of the model in the presence of local potentials will be proposed 
in this section, too. Methods of adding the potential to the soliton equation 
of motion will be described in section 3. Interaction of solitons with the 
potential walls and potential wells will be investigated in section 4. Some 
conclusions and remarks will be presented in the final section. 

2 Nonlinear Klein-Gordon Equation 

Different types of nonlinear Klein-Gordon equation have been proposed. 
They are different in their nonlinear terms which arise from application 
manners. Some of them are as follows: 
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- /3cP^ + jcP'' 



(2.1) 
= 0, (2.2) 

in which a and (3 are arbitrary constants and 7 = j^. Several locahzed 
solutions such as solitons, compactons and solitary waves have been found 
for these equations using different methods [U EJ [9l [10] . In this paper we 
will focus on the soliton solutions of the equation (2.2). The results are valid 
for the other equations, too. 

Wazwaz has found several localized solutions for the equation (2.2) using 
"Tanh" method. One soliton solution for this equation is [7] 
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where xq is the soliton initial position and u is its velocity. The signs "+" 



and "-" in the equation (2.3) denote to the kink and antikink solutions 
respectively. Figure 1 presents kink and antikink solutions with u = 0.5, 
a = 1, Xo = and Q; = /5 = latt = 0. 

As mentioned before, the Klein-Gordon type field theories have been 
applied to many investigations. The quantum aspects of the black hole 
have been studied considering a time-dependent classical black hole solution 
induced by the Klein-Gordon soliton [6]. The Hawking radiation and the 
Bondi reflected energy of the back reaction of the metric were calculated by 
adding a lump of energy through the Klein-Gordon field as 

f/(^) = f(^^-f)'- (2.4) 

The parameters j3 and fx are constant if soliton is sufficiently narrow. In 
this situation the effects of the gravity on the soliton characters are ne- 
glected. For broader energy lumps, we cannot take constant parameters for 
the soliton and its characters which are spatial functions|lHll2j. Therefore, 
a better situation can be investigated if we can consider a potential for the 
Klein-Gordon part of the theory with space dependent parameters like j3{x) 
and /u(x). 

The behaviour of transition regions in one-dimensional multistable con- 
tinuous media is described by the nonlinear Klein-Gordon equation with 




Figure 1: Soliton solution of the NKG described by equation (2.3). 



dissipation (NKGD) [B] 

//(Ait + 70t - D^xx + P{(t>) = 0, 



(2.5) 



where -?(</>) is a function of the field (p. Generally, the parameters of the 
model are some functions of the medium characters like temperature, vis- 
cosity and so on; however, such a system has not been studied with constant 
parameters yet. It is clear that we need a suitable model containing variable 
(at least slowly varying) parameters. 

Much work has been devoted to the dynamics of the one dimensional 
(ID) chain of coupled double- well oscillators and its applications in con- 
densed matter |14[ [T5| I16j . This model as well as the Frenkel-Kontorova 
chain became a very convenient and popular tool for the theoretical studies 
of dynamical and statistical properties of a number of materials with strong 
quasi- ID anisotropy |17j . One of the interesting systems of this class is a 
chain of hydrogen bonds where the proton transfers in adjacent hydrogen 
bonds are correlated because of the cooperative nature of hydrogen bonding. 
The proton dynamics in an ice crystal can be described using the discrete 
nonlinear Klein-Gordon (NKG) equationjlSl [T9] 



co</'xx + ^o^(0) = 0. 



(2.6) 



It is clear that a more realistic model contains disorders and disloca- 
tions in the hydrogen chain, which can be modeled using space dependent 
parameters in the NKG field theory. 

3 Two Models for "the NKG Soliton-potential" 
System 

Model 1: Consider a scalar field with the Lagrangian 

£ = ^a^</>5^</) - U {<P) (3.7) 

and the following potential: 

U (cP) = A(x) Qa<A' - ^/304 + ^7/) • (3.8) 

The equation of motion for the field becomes 

5^,5'^ (/) + A(x) («(/)- /3(/)3 + 705) = 0. (3.9) 

The effects of the external potential can be added to the equation of motion 
by using a suitable definition for A(x), like A(x) = 1 + v{x) 1201 EH ES]- For 
a constant value of parameter {X{x) = 1) equation (|3.9|) reduces to (2.2) 



and therefore (2.3) is its exact solution. Thus (2.3) is used as an initial 



condition for solving (3.9) with a space dependent A(x) when the potential 
v{x) is small. 



Hamiltomian density of the Lagrangian (3.7) is 



-Hi = ^<A' + ^<A' + Xix) Qa</.2 _ 1^^4 ^ 1^^6^ (3 10) 

Model 2: The potential also can be added to the Lagrangian of the 
system, through the metric of background space-time. So the metric carries 
the medium characters. The general form of the action in an arbitrary 
metric is 

1= f c{cp,d^(P)V^rxdt, (3.11) 

where "g" is the determinant of the metric g^'^{x). A suitable metric in the 
presence of a smooth and slowly varying weak potential v{x) is 



S-M^l '^f °^ |. (3.12) 
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The equation of motion for the field </> which is described by the Lagrangian 



(3.7) in the action (3.11) is 
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(3.13) 



This equation of motion in the background space-time (3.12) becomes [26] 

d'^(t) 1 d'^cj) dU{(p) 



{l + v{x)) 



dt"^ 1 + v{x) dx"^ 



+ 



d<p 



0. 



The field energy density is 



-H2 = {g'\x))'\4>' + \4>'+g'\x)U{4>). 



(3.14) 



(3.15) 



The energy density is calculated by varying both the field and the metric 
(See page 643 equation (11.81) of |25].). 



Solution (2.3) can be used as initial condition for solving (3.14) when 



the potential v{x) is small. 



4 Numerical Simulations 

Several simulations with different types of external potential v{x) have 
been performed using two presented models. The smooth and slowly varying 
potential v{x) = ae~^^^~'^> has been used in simulations which are reported 
below. Parameter "a" controls the strength of the potential, "b" represents 
its width and "c" adjusts the center of the potential. If a > 0, the potential 
shows a barrier and for a < we have a potential well. 

Simulations have been performed using 4th order Runge-Kutta method 
for time derivatives. Space derivatives were expanded using finite difference 
method. Grid spacing has been taken h=0.01, 0.02 and sometimes h=0.001. 
Time steps have been chosen as \ of the space step "h" because of nu- 
merical stability considerations. Simulations have been set up with fixed 
boundary conditions and solitons have been kept far from the boundaries. 
We have controlled the results of simulations by checking the total energy 
as a conserved quantity during the simulation. 

A moving soliton has kinetic and potential energies. Consider a static 
(zero velocity) soliton, in this situation the soliton has only potential energy. 
The difference between the energies of the two static solitons in different po- 
sitions comes from the difference between the potentials in those places. 



This means that one can find the potential as a function of collective co- 
ordinate X using the energy of a static soliton in a different position X. 
The shape of the potential (as a function of X) has been found by placing 
a static soliton at different positions and calculating its energy. It is clear 
that its energy is equal to the potential in that position. Figure 2 shows the 
shape of the potential barrier v{x) = 0.5e~^^ in different models which has 
been calculated using this method. This figure shows that the static parts 
of the potential is almost the same in both models. But we will find that the 
dynamical behaviour of the solitons are very different in the above models. 
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Figure 2: Potential barrier v{x) 
models. 



0.5e as seen by the soliton in two 



Suppose a soliton is placed far away from the potential. It goes toward 
the barrier in order to interact with the potential. There exist two different 
kinds of trajectories for the soliton during the interaction with the barrier 
depending on its initial velocity, which are separated by a critical velocity 
Uc- In low velocities Ui < Uc, soliton refiects back and reaches its initial place 
with final velocity of uj ^ Ui. Figure 3 presents trajectories of a soliton with 
different values of initial velocity for model 2. Critical velocity can be found 
by sending a soliton toward the potential with different initial velocities. 

Figure 4 presents the critical velocity as a function of the potential height 
in two models. This figure shows that the general behaviour of the soliton 
in the two models is the same. However, there is a little difference between 
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Figure 3: Soliton trajectory during the interaction with the potential barrier 



v{x) = 0.5e' 



-4x^ 



with different initial velocities for model 2. 



the results of model 1 and 2. This originates from the difference between the 
rest mass of the soliton in two models. Let us compare the soliton energy 
in model 1 (equation (3.10)) and model 2 (equation (3.15)). Static parts 
of the energy in two models are almost the same while the kinetic energy 
in two models is different. The effective mass of the soliton in model 1 is 
Mgffi = 1 while the effective mass in model 2 is Mf,ff2 = (1 + v{x))'^. The 
soliton rest mass of model 2 is greater than the soliton rest mass of model 
1. The difference between the soliton rest mass in the two models increases 
when the potential height becomes greater. Thus the difference between the 
soliton critical velocity in two models become larger, which is clearly shown 
in figure 4. 

Simulations show that scattering of a soliton with a potential barrier is 
nearly elastic. The soliton radiates a small amount of energy during the 
interaction. The radiated energy during the interaction in models 1 and 2 
is almost the same. 

Scattering of topological solitons on a potential well is more interesting. 
Unlike a classical point particle which always passes through a potential well, 
a soliton may be trapped in a potential well with specific depth[26l|27j. It is 
shown that the (f)^ soliton has not a fixed mass during the interaction with 
a potential well |28j and it works for the NKG soliton too. So we cannot 
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Figure 4: Critical velocity as a function of the barrier height of the potential 

v{x) 



0.5e ^^ in two models. 



look at the soliton as a point particle in some cases. Several simulations 
have been done for the NKG soliton-well system using the two models. Like 
the case of potential barrier, the general behaviour of the system is almost 
the same in both models. However, there are some differences in details of 
the interactions. Figure 5 presents a comparison between the shapes of the 
potential well v{x) = — 0.5e~^^ as seen by a NKG soliton in two models. 
There is not a critical velocity for a soliton-well system, but we can define an 
escape velocity. Consider a soliton which is located in the initial position xq 
near the center of the potential v{x). As figure 5 clearly shows, the soliton 
energy is a function of its initial position. This soliton can escape to infinity 
if its initial speed becomes greater than an escape velocity. The escape 
velocity for this situation is defined as the one for which the asymptotic 
speed is null at infinity. A soliton with initial velocity above the escape 
velocity, passes through the well and a soliton with initial velocity lower 
than the escape velocity falls into the well and is trapped by the potential. 
The rest mass can be calculated with the integration of Homiltonian 



density (Eq.(3.10) for model 1, Eq.(3.15) for model 2) respect to the position 
X. The calculated rest mass using models 1 and 2 are very near to each 
other. The rest mass in two models are exactly the same because of the 
same static part of the energy in two models. Integration of Hamiltonian 




Figure 5: The shape of the potential as seen by the sohton for the potential 
v{x) = — 0.5e~^^ 

density is calculated by Rumberg method for different values of the potential 
height. 

As figure 5 shows the effective potential for these two models are the 
same. Therefore it is expected that the escape velocity of the soliton in 
models 1 and 2 finds the same values for different values of the potential 
height. Figure 6 presents the escape velocity of a soliton in the potential 
well as a function of the potential depth. However, the potentials in the 
models are the same, but the calculeted escape velocities using these models 
are not equal. The reason for this difference can be explained using figure 7 
as follows. 

Figure 7 presents the trajectory of a soliton with an initial velocity 
u=0.35 during the interaction with potential v{x) = — 0.5e~^^' . The fi- 
nal velocity of the soliton in model 2 is less than the final velocity of the 
soliton in model 1. This means that the energy loss due to radiation in 
model 2 is greater than model 1. This is the main reason for the differences 
between the critical velocities observed in figure 6. 

The most interesting behaviour of a soliton during the process of scat- 
tering on a potential well is seen in some very narrow windows of initial 
velocities. At some velocities less than the Uc the soliton may reflect back 
or pass over the potential while one would expect that the soliton should be 
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Figure 6: Velocity as a function of the depth of the potential well. 
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Figure 7: Trajectory of a soliton with initial velocity u=0.35 during the 
interaction with the potential well v[x) = — 0.5e~^^ in both models. 
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trapped in the potential well. These narrow windows can be found by scan- 
ning the soliton initial velocity with small steps. Figures 9 presents this phe- 
nomenon in models 1 and 2. Figure 8(a) shows that a soliton with an initial 
velocity within the window 0.062 < Ui < 0.065 reflects back after the inter- 
action with the potential v{x) = — 0.4e~^^ simulated using model 1. Figure 
8(b) presents the same situation simulated using model 2. We couldn't find 
a soliton reflection in model 2 for the potential v{x) = — 0.4e~^^ . Some 
simulations have been performed for different values of the potential height 
(a=-0. 2,-0. 3,-0. 5,-0. 6,-0. 7,-0. 8) in model 2, but such a soliton reflection from 
the potential well has not been observed at all. Model 2 is built by varying 
the Lagrangian density with respect to both " field" and " metric" , therefore 
an energy exchange between the field and the potential is possible in this 
model. Some authors have related such a soliton reflection to the energy 
exchange between the field and the potential \1S\ I29j . It seems that this 
point is not true for the NKG solitons. Our simulations for the NKG soli- 
tons did not confirm this situation for model 2. This phenomenon needs 
more investigations. 
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Figure 8: (a) Soliton reflection from the potential well v{x) 
model 1. (b) Soliton reflection from the potential well v{x) 
model 2. 
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5 Conclusion 

Two models have been used to study the interaction of the NKG solitons 
with defects. The general behaviour of a soliton during the interaction with 
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external potentials in two models are almost the same. Model 1 adds the 
potential to the equation of motion by a different method from the one 
that is used in the model 2. Model 1 is easy to simulation while model 2 
is more analytic. The two models confirm that the interaction of a soliton 
with a potential barrier is nearly elastic. At low velocities it reflects back 
but with a high velocity climbs the barrier and passes over the potential. 
Soliton radiates some amounts of energy during the interaction with the 
potential. There exists a critical velocity which separates these two kinds 
of trajectories. The interaction of a soliton with the potential well is more 
inelastic. It is possible that a soliton scatters on a potential well and reflects 
back from the potential. This phenomenon depends on the model. Model 1 
predicts this behaviour, but this is in contradiction with what can be seen 
for the Sine-Gordon and 0"^ models. These models show this behaviour for 
model 2 while such a behaviour is observed in model 1 for the NKG solitons. 
It is interesting to investigate scattering of solitons of other models on 
defects using model 2. These studies help us to gain more comprehensive 
knowledge about the general behaviour of solitons. 

References 

[1] M. Bawin, M. Jaminon, Pion form factor and the Klein-Gordon Equat- 
tion, Physical Review C 30 (1) (1984) 331-334. 

[2] Y.R. Kwon, F. Tabakin, Hadronic atoms in momentum space, Physical 
Review C 18 (2) (1978) 932-943. 

[3] F.S. Catahotti, S. Burger, C. Fort, P. Maddaloni, F. Minardi, A. Trom- 
bettoni, A. Smerzi, M. Inguscio, Josephson junction arrays with Bose- 
Einstein condensates. Science 293 (5531) (2001) 843-846. 

[4] F.S. Catahotti, L. Fallani, F. Ferlaino, C. Fort, P. Maddaloni, M. In- 
guscio, Superfluid current disruption in a chain of weakly coupled Bose- 
Einstein condensates. New Journal of Physics 5 (2001) 71.1-71.7. 

[5] Elgarayahi A. New periodic wave solutions for the shallow water equa- 
tions and the generalized Klein-Gordon equation. Commun Nonlin Sci 
Numer Simul 2008; 13(5) :877-88. 

[6] Sung-Won Kim and Won T. Kim, Phys. Lett. B361 (1995) 38-44, (g^ 
|qc/9506064] 

[7] Wazwaz A.M. (2006) Chaos, solitons and fractals 28 1005-1013 Com- 
pactons, solitons and periodic solutions for some forms of nonlinear 
Klein-Gordon equations 

13 



[8] Sirendaoreji. Exact travelling wave solutions for four forms of nonlinear 
Klein-Gordon equations Physics Letters A 363 (2007) 440-447 

[9] Mustafa Inc. New compacton and solitary pattern solutions of the non- 
linear modified dispersive Klein-Gordon equations Chaos, Solitons and 
Fractals 33 (2006) 1275-1284v 

[10] Elgarayhi A. New periodic wave solutions for the shallow water equa- 
tions and the generalized Klein-Gordon equation Communications in 
Nonlinear Science and Numerical Simulation 13 (2006) 877-888 

[11] Stuart D.M.A., J. Math. Pures Appl. 80, 1, (2001) 5183 

[12] Stuart D.M.A., Ann. Sci. cole Norm. Sup. (4) 37 (2004), 312362 

[13] Otwinowski M.,Tuszynski J. A., Dixon J.M. , Phys. Rev. A 45, 10 (1992) 

7263-7273 
[14] J. A. Krumhansl and J. R. Schrieffer, Phys. Rev. B 11, 3535 (1975). 
[15] Gupta N., Deb B.M, Chemical Physics 327 (2006) 351-360 
[16] Radzihovsky L., Gurarie V., Phys. Rev A 81 (2010) 063609 

[17] Coulomb J. P. ,Floquet N., Andre G., Microporous and Mezoporous 
materials 141 (2011) 43-48 

[18] Zolotaryuk A. V., Savin A. V., Economou E. N., Phys. Rev. B 57, 1 

(1998) 234-245 
[19] Gorbach A.V. and Johansson M. Phys. Rev.E 67 (2003) 066608 

[20] Kalbermann G. Decay to bound states of a soliton in a well Chaos, 
Sohtons and Fractals 12 (2001) 625-629 

[21] Javidan K. (2006) J. Phys. A: Math. Gen. 39 10565-10574 
[22] Al-Alawi J. and Zakrzewski W.J. (2007) J. Phys. A 40, 11319-1131 
[23] Kalbermann G. A model for soliton trapping in a well. Chaos, Solitons 
and fractals 12 (2001) 2381-2385 

[24] Javidan K. (2008) Phys. Rev. E 78, 046607 

[25] Felsager B. (1998) Geomrtry, Particles, and fields (New York: Springer- 

Verlag) 
[26] Fei Z, Kivshar Y S and Vasquez L (1992) Phys. Rev. A46 5214 

[27] Piette B.M.A.G and Zakrzewski W.J. (2007) J. Phys. A: Math. Theor. 
40 No 2,329-346 

[28] Hakimi E, Javidan K. (2009) Phys. Rev. E 80, 016606 
[29] Y. S. Kivshar, Zhang Fei and L. Vazquez, Phys. Rev. Lett. 67, 
1177(1991). 



14 



